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ABSTRACT 

We present results from a Chandra observation of the NGC 346 star formation region, which 
contains numerous massive stars, and is related to N66, the largest Hll region of the SMC. In 
this first paper, we will focus on the characteristics of the main objects of the field. The NGC 346 
cluster itself shows only relatively faint X-ray emission (with L^ abs ~ 1.5 x 10 34 erg s“ l ), tightly 
correlated with the core of the cluster. In the field also lies HD 5980, a LBV star in a binary (or 
triple system) that is detected for the first time at X-ray energies. The star is X-ray bright, with 
an unabsorbed luminosity of ~ 1.7 x 10 34 erg s” 1 , but needs to be monitored further to 

investigate its X-ray variability over a complete orbital cycle. The high X-ray luminosity may be 
associated either with colliding winds in the binary system or with the 1994 eruption. HD 5980 
is surrounded by a region of diffuse X-ray emission, which may be a superimposed supernova 
remnant. 

Subject headings: (galaxies:) Magellanic Clouds-X-rays: individual (NGC 346, HD 5980) 


1. Introduction 

The Small Magellanic Cloud (SMC) is a an ir- 
regular galaxy at a distance of 59 kpc (Mathewson 


^nstitut d’Astrophysique et de Geophysique, Universite 
de Liege, Allee du 6 Aout 17, Bat. B5c, B 4000 - Li&ge 
(Belgium); naze@astro.ulg.ac.be 

2 School of Physics &: Astronomy, University of 

Birmingham, Edgbaston, Birmingham B15 2TT (UK); 
jmh@star.sr.bham.ac.uk, irs@star.sr.bham.ac.uk 

3 Universities Space Research Association, High En- 
ergy Astrophysics Science Archive Research Center, God- 
dard Space Flight Center, Greenbelt, MD 20771; corco- 
r an @ bar n egat . gsfc . n as a . go v 

4 Astronomy Department, University of Illinois, 1002 W. 
Green Street, Urbana, IL 61801; chu@astro.uiuc.edu 

5 Instituto de Astronomia, Universidad Nacional 
Autbnoma de Mexico, Apdo. Postal 70-264, 04510, Mexico 
D.F. (Mexico); gloria@astroscu.unam.mx 

6 Departement de physique, Universite de Montreal, 
C.P. 6128, Succ. Centre-Ville, Montreal, QC, H3C 3J7 
(Canada); moffat@astro.umontreal.ca 

7 Facultad de Ciencias Astronbmicas y Geofisicas, 
Universidad Nacional de la Plata, Paseo del 
Bosque S/N, B19000FWA La Plata (Argentina); 
virpi@fcaglp.fcaglp.unlp.edu.ar 
8 Research Fellow F.N.R.S. 


et al. 1986) that forms a pair with the Large Mag- 
ellanic Cloud (LMC). Both are satellites of our 
own Galaxy. The interstellar extinction towards 
the Magellanic Clouds is low, allowing studies 
of the X-ray sources to be undertaken with only 
small contamination. 

Previous X-ray observations of the SMC have 
been made with the Einstein Observatory , A SC A 
and ROSAT. These observations included surveys 
of the point source population (Kahabka et al. 
1999; Schmidtke et al. 1999; Haberl et al. 2000; 
Sasaki, Haberl, h Pietsch 2000; Yokogawa et al. 
2000) and in particular studies of the properties 
of the X-ray binary population (Haberl & Sasaki 
2000). The recent launch of Chandra , however, 
provides the opportunity to study the SMC with 
a far greater sensitivity and spatial resolution than 
ever before. 

For this investigation of the SMC, the N66 
nebula was chosen (Henize 1956). N66 is a giant 
H II region excited by the NGC 346 cluster, which 
is the largest star formation area in the SMC. 
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This cluster contains numerous massive stars, 
of spectral type as early as 02 (Walborn et al. 
2002). But the large amount of massive stars is 
not the only feature of interest in this field. On 
the outskirts of NGC 346 lies the remarkable star 
HD 5980, which underwent a LB V-type eruption 
in 1994. This massive binary (or triple?) system 
was subsequently monitored and analysed for its 
varying spectral and photometric properties, in 
visible and UV wavelengths. However, HD 5980 
has never been detected by past X-ray observato- 
ries. Previous studies of the SMC have discovered 
a bright extended source along the same line-of- 
sight as HD 5980 that has been attributed to a 
SNR. For the SMC, the NGC 346 field thus pro- 
vides one of the best opportunities to conduct an 
investigation of the X-ray domain. 

In this paper, we will describe in § 2 the ob- 
servations used in this study, then analyse the 
available data on NGC 346, HD 5980, and the ex- 
tended emission in § 3, 4, and 5, respectively. 
Finally, we will conclude in § 6. The second part 
of this work (paper II) will describe the properties 
of the other sources detected in the field. 


2. The Observations and Data Analysis 
2.1. X-ray observations 

NGC 346 was observed with Chandra for the 
XMEGA consortium on 2001 May 15-16 for 
100 ks (98.671 ks effective, ObsID = 1881, 
JD~2 452 045. 2d). The data were recorded by 
the ACIS-I0, II, 12, 13, S2 and S3 CCD chips 
maintained at a temperature of -*120°C. The data 
were taken with a frame time of 3.241s in the faint 
mode. The exposure was centered on the cluster, 
with HD 5980 lying 1.9' to the north-east of the 
aimpoint on ACIS-I3. Our faintest source possess 
a flux of about 2 x 10 32 erg s -1 (see paper II), 
assuming a distance of 59 kpc for the SMC. 

We excluded bad pixels from the analysis, us- 
ing the customary bad pixels file provided by the 
Chandra X-ray Center (CXC) for this particular 
observation. We have searched the data for back- 
ground flares, which are known to affect Chandra 
data, by determining the lightcurve of the to- 
tal count rate, but no flares were found. Event 


PI values and photon energies were determined 
using the FITS Embedded File (FEF) acisD2000- 
01-29fef_piN0001.fits. Note that this calibration 
matrix is still a preliminary product. 

For long exposures, removing the afterglow 
events can adversely affect the science analysis 
(underestimation of the fluxes, alteration of the 
spectra and so on) 9 . We thus computed a new 
level 2 events file by filtering the level 1 file and 
keeping the events with ASCA grades of 0, 2, 3, 
4 and 6, but without applying a status=0 filter. 
Throughout this paper, we will use this new file 
for all scientific analysis, except for the source de- 
tection, where it is better to use the pipeline level 
2 product to avoid mistaking afterglow events for 
real sources. 

Further analysis was performed using the CIAO 
v2.1.2 software provided by the CXC and also with 
the FTOOLS tasks. The spectra were analysed 
and fitted within XSPEC vll.0.1 (Arnaud 1996). 

2.2. HST WFPC2 images 

Two HST WFPC2 observations of the N66 
/ NGC 346 region in Ha are available from the 
STSci archives. They were taken on 1998 Sept. 27 
and 2000 August 7 for the programs 6540 and 
8196, respectively. HD 5980 is centered on the PC 
in the observations from program 6540, while pro- 
gram 8196 shows an area situated further south. 
The observations were made through the F656/V 
filter for 2070s (Prog. # 6540) and 2400s (Prog. 

# 8196). This filter, centered at 6563.7 A with 
a FWHM of 2 1.4 A, includes the Ha line and the 
neighboring [N II] lines, but we will neglect here 
the [N II] contamination. 

The calibrated WFPC2 images were produced 
by the standard HST pipeline. We process them 
further with IRAF and STSDAS routines. The 
images of each program were combined to remove 
cosmic rays and to produce a total-exposure map. 
Apart from HD 5980 (see Fig. 1), no obvious asso- 
ciation between X-rays and the Ha emission has 
been detected in the HST data. Fig. 1 will be 


9 see caveat on http://asc.harvard.edu/ciao/caveats/acis_cray.html 
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discussed more extensively in § 5. 

3. The NGC 346 Star Cluster 

NGC346 is the largest star formation region of 
the SMC. It contains a large fraction of the early 
O-type stars of this galaxy (Massey, Parker & 
Garmany 1989). It is related to the giant Hu re- 
gion N66 (Henize 1956), the largest and the most 
luminous one of the SMC (Ye, Turtle, & Kenni- 
cutt 1991). From many points of view, this ob- 
ject constitutes a unique feature in the SMC. The 
X-ray emission from the cluster itself has never 
been detected by previous X-ray observatories. 
Even though part of NGC 346 lies in the ACIS-I 
CCD gaps, diffuse emission associated with the 
cluster is seen in the upper corner of ACIS-I3. 
This position corresponds to the core of the clus- 
ter, which contains very early massive stars. In 
this area, the source detection algorithms found 
only one source, coincident to the position of star 
MPG435 (Massey, Parker & Garmany 1989) or 
N26 (Niemela, Marraco, & Cabanne 1986). How- 
ever, the emission is not restricted to this star. 
Most of it correlates well with other bright blue 
stars, i.e. MPG 368, 396, 417, 470 and 476 (see 
Fig. 2). But as the X-ray emission of the indi- 
vidual components can not be easily disentangled, 
we decide to analyse only the overall emission of 
these sources from now on. 

The total count rate of this emission in the 
0.3 — 10.0 keV band is 3.04 ± 0.36 x 10 -3 cnts s~ 1 . 
If we define 5, M and H as the count rates in 
the 0.3 - 1.0 keV, 1.0 - 2.0 keV and 2 - 10.0 keV 
energy bands, as in Sarazin, Irwin, & Bergman 
(2001) and Blanton, Sarazin, & Irwin (2001), 
we can compute two hardness ratios (HRs), H 21 
and i/31, defined as ( M — S)/(M 4- S) and 
(H — S)/(H + S), respectively. The hardness 
ratios H 21 and i/31 of the cluster are 0.41 ±0.11 
and —0.01 ± 0.24, respectively. These ratios are 
comparable to the softer group of all detected 
point sources (see paper II). 

The spectrum of the NGC 346 cluster was ex- 
tracted using a carefully chosen background win- 
dow, as close as possible to the cluster. As 
we regard the cluster X-ray emission to be ex- 


tended, we have used the calcrmf and calcar} 
tools to generate the appropriate response files. 
The best fit (xl =l 0.90, iV^ 0 /=5 0, Z — 0.1 Z 0 ) to 
this spectrum was an absorbed mekal model with 
the following properties (see Fig. 3): N(H) — 
0.42q'i 7 x 10 22 cm -2 and kT = I.OIqJJ keV. 
The luminosity in the 0.3 — 10.0 keV band was 
5.5 x 10 33 erg s -1 , which corresponds to an un- 
absorbed luminosity of 1.5 x 10 34 erg s l . The 
hydrogen column is comparable to the value ex- 
pected for NGC 346 (see paper II), but the tem- 
perature is one of the lowest we found, except for 
the extended source. 

In the region of diffuse emission, 30 blue stars 
have been detected by Massey, Parker & Gar- 
many (1989): MPG 435, 342, 470, 368, 476, 396, 
487, 417, 370, 471, 467, 495, 451, 330, 455, 454, 
500, 445, 468, 375, 508, 395, 439, 371, 485, 561, 
374, 366, 557, 400 (by decreasing magnitude). Of 
these, 16 have known spectral types that we can 
use to convert magnitudes to bolometric lumi- 
nosities. The bolometric correction factors were 
taken from Humphreys & McElroy (1984). We 
can then predict the expected X-ray luminosi- 
ties using log{L^p abs ) = 1.131og (Lbol) — 11.89 
(Berghofer et al. 1997). If the X-rays were com- 
ing only from these 16 stars, we should expect a 
total unabsorbed luminosity of 2 x 10 33 erg s -1 . 
This represents only 13% of the detected lumi- 
nosity. Several hypotheses could explain this dis- 
crepancy: a metallicity effect (Berghofer’s relation 
was determined for Galactic stars), X-ray emission 
from the other (unclassified) stars in this region, 
additional emission linked to interactions in bi- 
nary systems (such as colliding winds), emission 
from wind-blown bubble(s) around star(s) of the 
cluster core, or an extended region of hot gas as 
the winds from the individual stars combine to 
form a cluster wind (Canto, Raga, &c Rodriguez 
2000). 

Chandra has also observed X-ray emission from 
other very young (age < 5 Myr) stellar clusters, 
comparable to NGC 346. The Arches clusters in 
the Galactic Center region has been observed to 
have emission from several different regions, with a 
total X-ray luminosity of L£ no6s ~ 4 x 10 35 erg 
(Yusef-Zadeh et al. 2002). NGC 3603 has also 
been observed with Chandra (Moffat et al. 2002) 
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and the total unabsorbed cluster luminosity is 
j^unabs ^ | x 10 35 erg s _1 (though in this case it 
is clear that a substantial fraction of the emission 
is from point sources in the cluster). For both 
the Arches and NGC3603 clusters the fitted tem- 
perature of the X-ray emission is higher than the 
kT ~ 1 keV value for NGC346. 

It is worth noting that both the NGC3603 
and Arches clusters are probably a little more 
massive than NGC 346 (comparing the estimated 
number of ionizing photons quoted for NGC 346 
by Massey, Parker & Garmany (1989) and for 
NGC 3603 and the Arches by Figer, McLean , 
& Morris (1999)), and presumably the total 
mass and energy injection by the stars in the 
NGC 346 cluster is correspondingly only some- 
what smaller. NGC 346 is rather less luminous 
than both NGC 3603 and the Arches. This degree 
of variation is somewhat surprising as, on theoret- 
ical grounds, we would expect the cluster X-ray 
luminosity to scale as , with M cl the integrated 
mass-loss rate from the stars in the cluster (Canto, 
Raga, & Rodriguez 2000), which we might expect 
to be broadly similar. It is perhaps worth not- 
ing that in the case of NGC 346, because of the 
detector orientation and the lower flux levels, it 
could be that we are underestimating the cluster 
luminosity. 

4. The LBV HD 5980 

If NGC 346 contains a large fraction of the mas- 
sive stars of the SMC, at least one of the stars in 
this field, HD 5980, is certainly unique in itself and 
worth further examination. This eclipsing binary 
(or triple system?) was classified as WN-hOB 
before 1980, but its spectral type changed to 
WN3+WN4 in 1980-1983, and to WN6, with no 
trace of the companion, in 1992. In 1994, this star 
underwent a LBV-type eruption, and presented 
at the same time a WN11 type. The eruption 
has now settled down and the values are coming 
back to the pre-eruption ones (for more details 
see Koenigsberger et al. (2000) and references 
therein). 

Chandra is the first X-ray telescope to detect 
this peculiar system. However, even with a lOOks 
exposure, the data possess a rather low signal- 


to-noise ratio. HD 5980 has a count rate of only 
2.98 ± 0.19 x 10~ 3 cnts s” 1 in the 0.3 - 10.0 keV 
band and a hardness ratio (#21, #31) of (0.43, 
0.34), a rather low value compared to the other 
detected point sources (see paper II). 

We have extracted a spectrum of HD 5980 using 
the CIAO tool psextract. An annular background 
region around the source was chosen, in order 
to eliminate contamination from the surrounding 
extended emission (see § 5). This spectrum is 
shown in Fig. 4a. The best fit model to the data 
has the following parameters: N(H) = 0.22$ ^ x 
10 22 cm -2 and kT = 7.04 4 3 3 3 5 keV for an absorbed 
mekal model and N(H) = 0.28$;Jf x 10 22 cm“ 2 
and r=1.74j;.g9 f° r an absorbed power-law. These 
hydrogen columns are consistent with the value 
expected for NGC 346 (see paper II). The de- 
rived absorbed X-ray luminosity of HD 5980 is 
Lx = 1.3 x 10 34 erg s -1 in the 0.3 - 10.0 keV 
band for both models. Using the normalisation 
factor of the mekal model, we found a volume 
emission measure of ~10 57 cm -3 for this star. 

We can try to compare HD 5980 with other WR 
stars detected in X-rays. The unabsorbed X-ray 
luminosity of HD 5980 is 1.7 x 10 34 erg s~* in the 
0.3 - 10.0 keV energy range and 9 x 10 33 erg s~ 2 in 
the ROSAT range, i.e. 0.2 - 2.4 keV. This places 
HD 5980 amongst the most X-ray bright single 
WN stars (Wessolowski 1996) and the brightest 
WR-bOB binaries (Pollock 2002) of the Galaxy. 

Two factors could explain this high X-ray lu- 
minosity. First, the fast wind from the post- 
eruptive phases (from 1300 kms" 1 in 1994 to 
~2000 kms -1 in 2000) should now collide with 
the slow wind (ejected with a velocity as low as 
200 kms"' 1 during the eruption), which will in- 
crease the post-eruptive X-ray luminosity. X-ray 
emission from colliding ejecta has already been 
observed for another LBV, i.e. 7) Carinae. Close 
to this stax, the detected X-ray emission is corre- 
lated to the high velocity ejecta of the Homunculus 
nebula (Weis, Duschl, & Bomans 2001), which 
was created after the last great eruption. For 
HD 5980, however, the ejecta from the 1994 erup- 
tion has not had sufficient time enough to form a 
detectable LBV nebula around the star, and no 
{ Homunculus-like’ nebula from a previous eruption 
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is visible around HD 5980 (see Fig. 1): any X-ray 
emission from the colliding ejecta should thus still 
be blended with the stellar emission from HD 5980 
and that may contribute to explain the high X-ray 
luminosity. Unfortunately, no X-ray observation 
of HD 5980 before the 1994 eruption is available, 
which would allow us to quantitatively compare 
the pre- and post-eruption luminosities and con- 
firm this predicted luminosity enhancement. 

The fact that HD 5980 is a close binary sys- 
tem, containing two very massive stars suggests 
that colliding winds may provide another source 
for the observed X-rays, and we can make an 
estimate of the expected X-ray luminosity in 
this case. For the system, we assume that 
the stellar winds of the stars have returned to 
their pre-eruption values and we assume that 
Ma= 1.4x 10 -5 M © yr“ 1 , u oo (i4)=2500 km s^ 1 
for the O-star, and Mq~ 2xl0“ 5 Mq yr -1 , 
v oo ( J B)= : 1700 km s" 1 for the WN Wolf-Rayet star 
(Moffat et al. 1998). Thus the momenta of each 
stars wind are very similar, and the total wind 
kinetic energy is 5 x 10 37 erg s— 1. These wind 
parameters lead to a momentum ratio of t] — 0.97 
(Usov 1992) and the wind-wind collision shock 
will lie halfway between the stars. In this config- 
uration we expect about 1/6 of the wind kinetic 
energy to pass perpendicularly through the two 
shocks and be thermalised (Stevens, Blondin, & 
Pollock 1992), and because the systems are rela- 
tively close much of this energy will be radiated 10 . 

Consequently, we would expect the X-ray lumi- 
nosity of HD 5980 to be ~ 10 36 erg s _1 or more, 
rather than the observed value of ~ 10 34 erg s~*. 
The discrepancy between the predicted and ob- 
served X-ray luminosities could be due to a range 
of factors, e.g. the winds may not be at the pre- 
eruption levels - lower mass-loss rates or lower 
wind velocities will translate to a lower X-ray lu- 
minosity. Alternatively, the wind momenta may 
not be nearly equal, as a low value of r/ also tends 
to lower the luminosity. 

Using the ephemeris of Sterken & Breysacher 


(1997), we have computed a phase of 4> = 

0.24 — 0.30 for our Chandra observation: this is 
close to the eclipse of star A (which is the 1994 
eruptor), by star B (</>=0.36). If the eccentric- 
ity is e ~ 0.3 a simple adiabatic model predicts 
that the X-ray luminosity will vary inversely with 
separation and so we might expect a change in 
the intrinsic X-ray luminosity of ~ 2 through the 
orbit and of ~ 10% during the Chandra observa- 
tion. We have detected a variation of apparently 
larger amplitude in the count rate of HD 5980 (see 
Fig. 4b). This still suggests we are seeing orbital 
variability associated with colliding winds, since 
the additional effects of changing absorption on 
the observed luminosity can either reduce or en- 
hance this variability. This observed variability is 
much more likely to be due to colliding winds than 
wind-blown bubble type emission, where no short 
timescale variability would be expected. The high 
fitted temperature of HD 5980 kT ~ 7keV also 
suggests that we are likely seeing colliding wind 
emission, as this temperature corresponds to a 
pre-shock velocity of ~ 2500 km s -1 , comparable 
to that anticipated. 

However, the detailed characteristics of the X- 
ray properties of HD 5980 need to be studied fur- 
ther, in conjunction with detailed colliding wind 
models (c.f. the case of T] Carinae; Pittard & 
Corcoran 2002). An XMM-Newton observation 
is scheduled at phases <f> = 0.09 - 0.10, just after 
periastron: the comparison between both datasets 
may enable us to detect further the signature of 
the colliding wing region, and constrain it more 
precisely. 

Finally, we may note that the unabsorbed X- 
ray luminosity of the central source in r) Carinae, 
if we assume a distance of 2.3 kpc, varies between 
0.6 and 2.5 xlO 35 erg s -1 (Ishibashi et al. 1999), 
which is much lower than HD 5980. 

5. The Extended Emission around HD 5980 

Around the position of HD 5980 lies a region of 
X-ray bright extended emission. It was first de- 
tected by the Einstein Observatory (source IKT 18 
in Inoue, Koyama, & Tanaka (1983)), and subse- 
quently observed by ROSAT (source [HFP2000] 


10 In terms of the cooling parameter \ defined by Stevens, 
Blondin, &; Pollock ( 1992 ), for the WN wind x < 1 and 
for the O-star wind \ ~ F 
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148 11 in Haberl et al. (2000)). The overall 
shape of the emission is more or less rectangular, 
with an extension to the northeast. It contains 
a few bright or dark arcs, but apart from these, 
the brightness is rather uniform, with no obvious 
limb-brightening. It has a size of 130" x 100", i.e. 
37pcx29pc at the SMC’s distance. HD 5980 lies 
towards the top center of the source. An X-ray 
bright filament extends from 9" west of the star to 
23" to its south. An X-ray dark feature appears 
some 30" at the southeast of HD 5980. The total 
count rate of this source in the 0.3 - 10.0 keV 
energy range is 7.25 ± 0.17 x 10~ 2 cnts s _1 , and 
the hardness ratios H21 and 7731 are -0.24 ±0.02 
and -0.96 ± 0.05, respectively: it is the only soft 
source present in the field. 

The spectrum of the extended emission was 
extracted with a carefully chosen background re- 
gion, located as close as possible to the source 
and in a region on the ACIS-I1 CCD where there 
are no point sources. As the X-ray emission is 
extended the rmf and arf matrices were calculated 
using the calcrmf and calcarf tools. The best fit 
(X^ = l*03, Nd 0 f= 244) to this spectrum was an ab- 
sorbed mekal model with the following properties 
(see Fig. 5): N(H) = 0.12°;}* x 10 22 cm' 2 , kT = 
O-66 q ;04 keV and Z = 0.17q The flux in the 

0.3 — 10.0 keV band was 3.35 x 10 -13 erg cm~ 2 
s“\ i.e. a luminosity of 1.40 x 10 35 erg s” 1 . The 
low value of the hydrogen column suggests that 
this X-ray source lies between NGC346 and the 
observer, but still in the SMC (see paper II). The 
hydrogen column is, however, larger than the value 
predicted using ultraviolet absorptions (Koenigs- 
berger et al. 2001; Hoopes et al. 2001). 

Using the normalisation factor of the mekal 
model, we found a volume emission measure 
f n e Ti}{dV of ~3xl0 58 cm“ 3 . Considering a con- 
stant density for the hot gas, a sperical geometry 
of diameter ~33 pc for the extended emission and 
assuming a pure H composition, the density of the 
hot gas is roughly 0.2 cm -3 and the total mass of 
the hot gas is ~100 M 0 . 

Using the available data, a deeper analysis was 


11 This ROSAT source may encompass some of the X-ray 
emission from the cluster. 


conducted, searching for the existence of a temper- 
ature gradient throughout the source area, with 
color and temperature maps. For color maps, 
we generated a set of narrow-band images, e.g. 
0.3 - 0.9 keV, 0.9- 1.2 keV and 1.2 - 2.0 keV, that 
enables us to distinguish the harder components 
from the softer parts. In addition, temperature 
maps were constructed from spectral fits in small 
regions of the SNR. The only obvious trend in both 
the color and temperature maps is the softness of 
the northeastern extension of the SNR. However, 
it is not possible to draw any firm conclusions re- 
garding the presence of any temperature gradient. 

We have also constructed images of the ex- 
tended emission in narrow wavelength bands, each 
corresponding to a specific ion. The data in each 
band were first binned to obtain 4.9" x 4.9" pixels, 
and then smoothed by convolution with a Gaus- 
sian of a = 0.5 pixel (see Fig. 6). In contrast to 
N132D (Behar et al. 2001), no significative differ- 
ences between the morphology of highly ionized 
species (Mg 10 + Ne 9 +, Si 12 +, Fe 19+ ) and lower 
ionization species were detected. 

5.1. The Nature of the Extended Emission 

Since its discovery the extended X-ray emis- 
sion has been attributed to a supernova remnant 
(SNR). Further evidence was subsequently ob- 
tained to support this hypothesis; for example, 
Ye, Turtle, & Kennicutt (1991) found a non- 
thermal radio source at this position thay they 
called SNR 0057 — 7226. Using the lowest con- 
tours of Ye, Turtle, & Kennicutt (1991), the 
radio source is 56 pcx52 pc, slightly larger than 
the X-ray source. 

Moreover, evidence of high velocity motions in 
this region were observed in the visible and UV 
ranges. Using echelle spectra centered on Ho, 
Chu & Kennicutt (1988) found clumps moving at 
+170kms -1 from the main quiescent component 
(vlsr ~ 300 kms -1 ): only the receding side of the 
expanding object is seen, suggesting a low density 
where the approaching side arises. UV analyses 
with IUE (de Boer & Savage 1980), HST STIS 
(Koenigsberger et al. 2001) and FUSE (Hoopes 
et al. 2001) have also confirmed the presence 
of an expanding structure. But while de Boer 
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& Savage (1980) and Hoopes et al. (2001) 
found only absorptions at V£SR =+300km s _1 , 
Koenigsberger et al. (2001) detected both ex- 
panding sides of the object, with components at 
vlsr = +21, -(-52, +300, +331 km s -1 . However, 
the absorptions at +21 and +52 km s'" 1 were 
small, and needed a future confirmation (Koenigs- 
berger, private communication). Absorptions at 
any of these velocities are not seen in the spectra 
of Sk 80, a close neighbor of HD 5980 situated on 
the edge of the X-ray source. A fast expanding 
structure close to the position of HD 5980 is thus 
present. 

Unfortunately, even the high-resolution HST 
WFPC2 images of the close neighborhood of 
HD 5980 do not show any clear nebula associ- 
ated with this object (see Fig. 1). However, a lot 
of filamentary structures - generally indicative of 
a SNR (Chen et al. 2000) - are seen throughout 
the field, but are not limited to the exact loca- 
tion of the X-ray source. There is thus no obvious 
Ha feature directly associated with the extended 
emission. 

Considering all the evidence (non-thermal emis- 
sion, high velocity expanding shell etc), the X-ray 
/ radio source should thus be regarded as a SNR 
located in front of HD 5980 but still belonging 
to the SMC. The SNR hypothesis is further sup- 
ported by the comparison of the size of this fea- 
ture to its X-ray and radio fluxes (Mathewson et 
al. 1983). 

However, as HD 5980 is projected more or less 
in the middle of this X-ray emission, it is hard 
to believe that there is no link at all between the 
two and that we have here only a chance super- 
position, especially when we know that HD 5980 
has undergone a LBV-type eruption in 1994. This 
1994 eruption is unlikely to be the only one un- 
dergone by the system. Can we find any evidence 
showing that HD 5980 is directly responsible for 
this emission? 

First of all, we can compare HD 5980 with other 
objects, since there are known interactions of other 
LBVs with their surroundings, e.g. r) Carinae, 
the most famous LBV of the sky. This star is 
indeed surrounded by an X-ray bright extended 


emission associated with the entire Carina Neb- 
ula: it is comparable in X-ray temperature (kT), 
size and morphology to the X-ray emission around 
HD 5980 (Seward & Mitchell 1981 ; Fig. 7). But 
it is generally assumed that this Nebula has been 
formed through the collective action of all stars 
of the Carina cluster, not only by the single ac- 
tion of 7 ] Carinae. Also, the expansion velocity of 
the Carina Nebula is far smaller than ours: 15- 
20 kms -1 , see Walborn &; Hesser (1975). All 
high velocity structures in this nebula are smaller 
structures localized closely on a few stars. 

However, Walborn (1978) claimed to have 
found some evidence of interaction between HD 5980 
and its surroundings: an arc of radius ~ 30" cen- 
tered on the star is actually visible in Ha and 
[O ill] (but not in [Sii]). But this arc is rather 
diffuse (see Fig. 1), indicating a lack of shock 
compression. Moreover, its velocity is the one of 
the quiescent region, and no splitting region is 
seen beginning at the arc and extending towards 
the star: this arc is certainly not the rim of an 
expanding shell. Finally, no strong correlation 
with the X-ray bright arc can be derived: both 
have different sizes and shapes, as can be appreci- 
ated on Fig. 1, and we can therefore conclude that 
they probably belong to two distinct features. The 
comparison between X-rays and visible data thus 
do not show any clear indication of an interaction 
between HD 5980 and its environment. 

Other clues also seem to reject such a possi- 
bility. For an interaction between a LBV neb- 
ula and circumstellar material surrounding the 
star, we should detect an enhancement of the 
[N Ii]/Ha ratio, like for SNR 0540 — 69.3 in the 
LMC, SNR 0102 — 72.3 in the SMC or the ultra- 
luminous SNR in NGC6946 (Dunne, Gruendl, & 
Chu 2000): this is not seen in the echelle spec- 
tra of this area and we have therefore to reject 
this possibility. We also note that none of these 
interactions could explain the redshifted UV ab- 
sorptions in the spectra of HD 5980, nor the non- 
thermal radio emission from the extended source. 

The final possibility is that the source is in- 
deed a SNR, but that HD 5980 is still directly 
responsible for it. The third component of this 
system could be a compact object, a remnant of 
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the star that exploded as a SN long ago. However, 
as the SNR is now in front of HD 5980 (as indi- 
cated by the UV absorptions and the low hydrogen 
column), HD 5980 should have been ejected and 
have velocity of > +300 km s -1 . Unfortunately, 
Koenigsberger et al. (2000) show that the aver- 
age velocity of the system is entirely compatible 
with that of the NGC 346 cluster (but this conclu- 
sion only applies to the radial component, since 
nothing is known on the tangential motion of the 
star). Therefore we conclude that, most probably, 
another very massive star from the NGC 346 field 
- whose stellar remnant has never been detected 
up to now - is the probable progenitor of the SN 
explosion. 


6. Summary 

In this serie of articles, we report the analysis 
of the Chandra data of NGC 346, the largest star 
formation region of the SMC. In this first paper, 
we have focused on the most important objects of 
the field : the NGC 346 cluster and HD 5980. 

The cluster itself was relatively faint, with a 
total luminosity of L^p abs ~ 1.5 x 10 34 erg s“*. 
Most of this emission seems correlated with the 
location of the brightest stars of the core of the 
cluster, but the level of X-ray emission probably 
cannot be explained solely by the emission from 
individual stars. 

In this field lies another object of interest: 
HD 5980, a remarkable star that underwent a 
LBV- type eruption in 1994. Chandra is in fact 
the first X-ray telescope to detect HD 5980 indi- 
vidually. In X-rays, the star appears very bright, 
comparable only to the brightest WR stars in the 
Galaxy. The comparison of our results with future 
X-ray observations will enable us to better under- 
stand HD 5980: for example, phase-lock variations 
would be analysed in the perspective of the col- 
liding wind behaviour of this binary, while other 
variations may be related to the recent LBV erup- 
tion. Follow-up observations are thus needed to 
complete the study of this system. 

A bright, extended X-ray emission is seen to 
surround this star, but it is probably due to a SNR 


unrelated to the star, the progenitor of which is 
unknown. We note however its close resemblance 
to the Carina Nebula in which the LBV rj Carinae 
lies. 
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Fig. 1. — HST WFPC2 Ha data of the region 
surrounding HD 5980 (which is the star located at 
the center of the PC chip), with the adaptively 
smoothed Chandra X-ray data overlaid as con- 
tours. The extended X-ray emission and the X- 
ray emission associated with HD 5980 are clearly 
seen. Parts of the NGC346 cluster are visible in 
the SW portion of the field. 


Fig. 3.—- The X-ray spectrum of the NGC346 
cluster (Src 36), shown along with the best-fit 
absorbed mekal model, with N(H) = 0.42 x 
10 22 cm -2 and kT = 1.01 keV (see § 3). 
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HD 5980 + NGC346 
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Fig. 2. — The smoothed Chandra X-ray data of the 
NGC 346/HD 5980 region (smoothing by convolu- 
tion with a gaussian of a = 2.5") superimposed 
on the DSS image, showing the extended nature 
of the X-ray source associated with the cluster. 
The location of the brightest stars in the cluster 
are indicated. 



Fig. 4. — a. The Chandra X-ray spectrum of 
HD 5980, shown along with the best-fit absorbed 
mekal model, with N(H) = 0.22 x 10 22 cm” 2 and 
kT = 7.04 keV. b. The Chandra X-ray lightcurve 
of HD 5980, in the 0.3-10 keV range and with 10 
bins of 10 ks each. 
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SNR 0057-7226 



Fig. 5. — The Chandra X-ray spectrum of the 
extended emission around HD 5980, shown along 
with the best-fit absorbed mekal model, with 
N(H) = 0.12 x 10 22 cm” 2 , kT = 0.66 keV and 
Z = 0.17Z q (see § 5). 
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Fig. 7. — The Carina Nebula compared to the 
extended emission around HD 5980. If we as- 
sume a distance of 2.3 kpc for the Carina Neb- 
ula and 59 kpc for the SMC, both images are 
~55 pcx55 pc. ROSAT PSPC support rings and 
Chandra ACIS-I CCD gaps can be spotted on 
these images. Both images were smoothed by con- 
volution with a gaussian (cr=2.5" for HD 5980 and 
15" for the Carina Nebula). 


Fig. 6. — Narrow-band images of the extended 
emission surrounding HD 5980. Each image is la- 
beled with the principal line-emitting ion and the 
energy range used to generate the image. 
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